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Abstract. We present models for dusty high-redshift Lya emitting galaxies by com- 
bining the Press-Schechter formalism with a treatment of inhomogeneous dust distri- 
bution inside galaxies. These models reproduce the surface density of emitters inferred 
from recent observations, and also agree with previous non-detections. Although a 
detailed determination of the individual model parameters is precluded by uncertain- 
ties, we find that (i) the dust content of primordial galaxies builds up in no more than 
~ 5 x 10 8 yr, (ii) the galactic HII regions are inhomogeneous with a cloud covering 
factor of order unity and (iii) the overall star formation efficiency is at least ~ 5%. 
Future observations should be able to detect Lya galaxies upto redshifts of z ~ 8. If 
the universe is reionized at z T 55 8, the corresponding decline in the number of Lya 
emitters at z <; z t could prove to be a useful probe of the reionization epoch. 

INTRODUCTION 

Until a decade ago, the search for high-redshift Lya galaxies had enjoyed no 
compelling successes [1]. With the improved sensitivity on large area telescopes 
in the last couple of years, these young galaxies are finally being detected [2,3]. 
Clearly, this population is of great interest to the field of galaxy formation and the 
early evolution of the universe. Fundamental questions we need to understand are 
why the earlier surveys have been unsuccessful, how many objects are still expected 
to be found, and what physical conditions pertain in these early systems so that 
the Lya radiation may escape. 

In a galactic setting, the emerging luminosity of the Lya line is strongly modu- 
lated by the amount and spatial distribution of stellar dust [4]. Observations have 
firmly established the strong decrease in Lya equivalent width with increasing oxy- 
gen abundance [5]. Indeed, even a modest amount of dust (observationally traced 
by oxygen) inside a homogeneous HII region is sufficient to attenuate all of the pro- 
duced Lya radiation, because Ly alpha photons are resonantly scattered, and dust 
absorption significantly increases the effective line optical depths [6]. However, this 
situation can be alleviated if the medium is inhomogeneous [4]. In a multi-phase 
medium where dust resides in opaque clumps, the photons do not penetrate the 
clumps and spend most of their time in the interclump medium, where the opac- 
ity is very small. The result is that the escape of Lya radiation is significantly 
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FIGURE 1. The escape fraction F csc of Lya photons from an inhomogeneous medium, as a 
function of dust content. The three curves correspond to three different values of the covering 
factor of opaque clumps, F cov =1, 5, and oo. 

enhanced compared to the case of resonant scattering in a homogeneous medium. 
Alternatively, large enough velocity gradients can provide a second way of efficient 
escape of Lya radiation, as suggested by recent observations of 4 nearby galaxies [7] . 

In a cosmo logical context, the formation of dust requires the presence of metals. 
Since the enrichment of the gas with metals is expected to build up only gradually 
during the star formation history of the universe [8], at high enough redshifts the 
enrichment level will be low. The correspondingly small dust content of the earliest 
galaxies would facilitate the escape of their Lya radiation. Such an interplay be- 
tween the dust enrichment of the IGM and the galactic environment could be the 
mechanism that leads to the observed surface density of high-z Lya emitters [9] . 



LYa EMISSION FROM INHOMOGENEOUS CLOUDS 

The radiative transfer problem which needs to be solved for the Lya line is 
well studied [10]. We have modeled individual galaxies with a range of masses 
for the ionizing stars, dust content, and inhomogeneity, using a numerical Monte 
Carlo approach [11]. We assumed a Scalo IMF for the spectral types of the central 
stars in the HII regions, ranging from 05 to Bl, with the stars distributed in a 
statistically homogeneous manner inside a percolating multi-phase medium. We 
adopted the formalism of Neufeld [4] for a multi-phase medium, and assumed a 
typical line width of Av ~ 8 km/s. We parameterized the multi-phase medium 
by opaque clumps embedded in an inter-clump medium of negligible opacity. The 
clump covering factor F cov , i.e. the average number of clumps along a line of sight, 
fixes the degree of inhomogeneity. The escape fractions F esc were computed on 
a grid of models with dust contents between = 10~ 2 — 1 solar, and covering 
factors F cov = l,5,oo. Figure 1 shows the escape fraction, and demonstrates its 
sensitivity to the covering factor. The inhomogeneous percolating slabs are much 
more transparent than the homogeneous ones; the difference around ~ lO^ 1 
solar is over an order of magnitude. 
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COSMOLOGICAL ABUNDANCE OF LYa EMITTERS 

In order to model the cosmological abundance of high-redshift Lya emitters, we 
assumed that the formation of dark matter halos follows the Press-Schechter [12] 
theory, and every halo with a virial temperature larger than 10 4 K forms a galaxy 
that goes through a Lya emitting phase [13]. A fraction of the gas was turned 
into stars, with a constant star formation rate (SFR) over a period of years. The 
SFR was related to the intrinsic Lya luminosity [14], assuming case B recombina- 
tion. We adopted the simplest assumption, i.e. that e* and t* both have the same 
constant values in each halo. The amount of dust produced and retained in each 
galaxy, Zd,iSM then increases linearly for t* years, after which it reaches the final 
value of Zd,iSM(£*) = 0.3 solar. Similarly, we assume that each galaxy deposits dust 
into the surrounding IGM at a constant rate for t+ years, and causes an enrich- 
ment of the surrounding regions within the intergalactic medium to Z(1,igm(z)- We 
normalize Zd,iGM to the cluster metallicity of ~ 0.3 solar at redshift z ~ 1. Note 
that our ^a,iGM is the average metallicity within the polluted regions, and not the 
universal average dust fraction of the IGM. Given a cosmology, the five parameters 
t+, F cov , e*, Z d> iQ M (z = 1), and ZdjsM^*) uniquely determine the number density of 
Lya emitters at any flux and redshift. Figure 2 shows the resulting surface density 
of emitters, as a function of redshift, in a flat ACDM model with a tilted power 
spectrum (tt ,tt A ,tt h ,h,a 8h -i,n) = (0.35, 0.65, 0.04, 0.65, 0.87, 0.96). We also indi- 
cate the observed surface density of emitters [2] at the two redshifts z = 3.4 and 
z = 4.5. 

Figure 2 has several interesting features. It demonstrates that the surface abun- 
dance strongly depends on the time over which the dust is produced (t star ), and 
on the ambient inhomogeneity of the HII regions that surround the ionizing OB 
stars (-F C ov)- Our models indicate that the dust content builds up in ^5x 10 s yr, 
the galactic HII regions are inhomogeneous with a cloud covering factor of order 
unity, and the overall star formation efficiency is at least ~ 5%. These numbers 
should be predicted by more complete and detailed models of galactic evolution, 
and will be useful discriminators between such models. The surface density of emit- 
ters is also a steep function of the detection threshold, a feature that makes our 
model consistent both with recent detections [2] and earlier upper limits [1]. Our 
models predict that the surface density changes relatively slowly with redshift to 
z 8, and Lya galaxies will be detectable, around the present flux threshold, upto 
redshifts as high as ~8. Recent detections of three new Lya emitters at z ~ 5.7 
support this conclusion [15]. However, if the universe is reionized at z re ion ^ 8, 
then the damping wing of Lya absorption from the neutral IGM would severely 
damp the Lya emission line [16]. This would render Lya emission from z ^ z rc i on 
undetectable, and cause a decline in the number of observed emitters beyond z TC i on . 
In addition, the neutral IGM would likely imprint a characteristic assymetry on 
the emission line profiles of emitters located near the reionization epoch. Provided 
this assymetry is measured, the disappearance of the Lya emitter population could 
be useful diagnostic signature of reionization. 
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FIGURE 2. The surface density of Lya emitters in our standard model (solid lines) with fluxes 
above different values of the detection threshold. The two data points are taken from [2]. For the 
fixed threshold Fq = 1.5 x 10~ 17 erg cm -1 s _1 7 the dashed lines show the surface density when the 
star-formation rate is increased or decreased by a factor of 10. Similarly, the dotted lines show 
the surface density when the covering factor is changed to F cov — 1, or oo. 
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